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Abstract 

We propose to calculate bosonic and fermionic determinants with some general 
field background, and the corresponding 1-loop effective actions by evaluating random 
walk worldline loops generated statistically on the lattice. This is illustrated by some 
numerical calculations for simple gauge field backgrounds and then discussed for the 
general case. 
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1 Introduction 



The discussion of bosonic and fermionic determinants in some background and of the 
related 1-loop effective action is most transparent in the worldline/relativistic particle 
formalism |IJ-@. The Euler-Heisenberg effective action of a constant electromagnetic field 
||], for instance, can be presented very elegantly in this way Q. Also for general gauge 
field backgrounds this formulation is very helpful allowing for an economical inverse mass 
(derivative) expansion to high order || . For background fields like an instanton, sphaleron 
or a bounce solution used in tunneling problems the evaluation of a derivative expansion || 
and alternatively the approach based on solutions for the eigenvalue spectrum Jjj are rather 
sophisticated. The nonabelian gauge field groundstate background in 3- and 4-dimensional 
QCD and in the hot electroweak theory is even less accessible to an analytical treatment. 
For instance, the discussion of a nonperturbative magnetic and tachyonic mass in 3- and 
4-dimensional SU(2) theory is very difficult even if one uses simplified assumptions about 
the gauge field vacuum M . Finally, we should notice that 2-loop effective actions are also 
accessible g |, fR]] 

In this paper we consider the numerical evaluation of quantum field theoretical de- 
terminants in the worldline formalism using lattice methods. We consider the bosonic 
determinant in a gauge field background 



T(A) = - log det(£>^ + m' z ) = -tr log(£> 2 + m 2 ) 

%r f d D x [ [Dx] tr P ex.pl - [ dr ( ^- + igA^ + m 2 

o T J Jx(o)=x(T)=x I Jo V 4 
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where the worldline sum (here already in the Euclidean formulation adequate for the lattice 
analysis) is over all closed paths and trp indicates path ordering. In proceeding to the 
lattice formulation of the problem it is natural to consider the loops on the same (usually 
hypercubic) lattice, on which the general gauge field configurations are generated. This 
appears to be difficult if one wants to enforce the exp{— J Q dr^-} weight for randomly 
chosen paths, as proposed in an interesting recent paper [[IljJ. However, as is well known, 
random walk paths automatically implement the above integral weight of the free particle 
propagation on a worldline p^] . The expression (|l]) can thus be discretized on a D- 
dimensional lattice Ad as 

T*«*(A)= £ C(A,x); C(A 1 x) = J2je- m2 ^(±) L £ tr ( ]J U t 

zeAj, L v J Wl{x)} \ieu, L ( x ) 

where dimensionfull quantities are understood as given in units of the lattice scale a 
according to their naive dimension and we have also discretized T using the r intervals 

T a 2 
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Here {u>l(x)} is the set of all closed lattice paths with length \lol\ = L obtained by 
random walk (RW) starting and ending at x, and we have written down explicitly the 
random walk "measure" [2D)~~ L (to be implemented by the actual RW procedure). The 
mass dependence e~ m T could be also obtained by a random walk weight but we prefer to 
keep this dependence explicit. XJ\ are link variables. 

On the lattice the expansion according to the loop lengths (§) is of course equivalent 
to the usual hopping parameter expansion for the logarithm of the determinant, both for 



bosons and for fermions (see, e.g., [12], 13], ]|14[|). This can be seen explicitly by rewriting 
(e.g., in the bosonic case) 



K L = {m 2 + 2D)- L ~ (2D)- L e~ m2 ^ (4) 
m 2 < 2D 

(and similarly for fermions). The number of different loops of length L increases exponen- 
tially with the length, therefore one usually must cut the hopping parameter expansion 
at a rather small L. We suggest here to use a statistical procedure for sampling loops in 
a wide interval of lengths, instead of summing over all possible loops up to a given order 
(see also (l]|). In our approach we use the random walk representation (||) to construct a 
(sub-) ensemble of loops of various lengths and evaluate the determinant with help of this 
sub-ensemble. By construction, the loops appear with the correct probabilities, and since 
this is a statistical procedure we may expect to have only statistical errors. Since loops 
of length L are typically relevant at the scale y/L, arbitrarily cutting the loop expansion 
at some limited L may completely miss physical effects at scales larger than x/~L. In our 
method we trade these systematic errors for statistical ones such that we may hope to 
reproduce also such effects, albeit within statistical uncertainties. Notice that, on the 
one hand, we cannot use too rough lattices since this introduces discretization errors. On 
the other hand, the contribution of large loops is damped even at small mass by a factor 
oc L~ D / 2 , see (|8|) next section. This suggests that, depending on the effects we want to 
observe, there is an optimal range of loops to be taken into account. In practice we use 
loops of length up to L ~ 200 on lattices of 32 D . For the cases studied here these ensembles 
appeared adequate. 

Chapter 2 describes the procedure how to generate closed loops and checks that the 
case without interaction on our loop ensembles converges to the well-known continuum 
(large L) result. In expression we can insert some (quasi) classical background gauge 
field An discretized by lattice connections U. This we will do for a constant magnetic 
field in chapter 3. This case allows comparison with the well-known Euler-Heisenberg 
action in its worldline formulation. We also treat there the case of constant SU{2) fields 
in two perpendicular 2-planes in four dimensions which has a topological meaning 
The case of a constant magnetic field in a half space treated in ref. [pTTfj , which introduces 
a correlation scale (penetration length) is also studied. 

For a quantum gauge field vacuum the action Y dlscr (A) has to be inserted into a gauge 
field path integral discretized on the lattice to obtain the (quenched) vacuum energy 

VG9) = <r> = £ \e- m ^ r_L) L y: ^ / [DU]e-f>8™m)tv ( n vA (5) 

l {ui L } ye^L / 
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with the Yang-Mills discretized action Sym an d partition function Zym- I n the case of a 
further field (ft in the background coupling to the loop a term V"((j>(x)) has to be introduced 
in the discretized worldline action and further weights the sum over random paths. These 
1-loop actions in a general background are our real goal to be accomplished in future work. 
Chapter 4 presents a final discussion and an outlook to the case of spin in the loop, in 
particular to fermionic determinants. 

2 The lattice loop ensemble 

As remarked in the introduction we use random walk to produce a loop ensemble ther- 
malized with the "Boltzmann factor" 

see (||), which is the discretized form of the corresponding factor in ([!]). In reconstructing 
the whole path integral we only need to weight the RW loop contributions tr (lliew L Ui) 
with 

L-V™ 2 ^ (7) 

since the factor (2D)~ L is implicit in the RW generation. From now on we denote by {lo}, 
{wi} the corresponding ensemble of loops as it is produced by the random walk procedure, 
that is, with the frequencies with which the loops are generated in the actual procedure. 

In generating the loops we specify a maximal number of trials M (e.g., 2000000 for 
the 32 3 lattice) for each given number of steps (between 4 and 180). We first generate 
an {uj(x^ = 1)} ensemble of loops by starting the random walk at the point x„ = 1 
of the periodic lattice and collecting the loops which after the prefixed number of steps 
have returned to x^ = 1, independently on possible self-crossings or retracing. We define 
£(A, 1) using this {w(x M = 1)} ensemble (for the application in sect. |3]^ we first shift 
the approximate "center of mass" of the loops to x^ = 1). For any given x we can 
calculate C(A, x) in (J2j) by accordingly shifting the loops of the {^(x^ = 1)} ensemble 
as a whole. On the other hand, the full loop ensemble, that is, the x-summation Y^xeA D 
in (g), is approximated by considering sufficiently many random translations of each of 
the loops of the {uj(x^ = 1)} ensemble (for homogeneous field configurations this is, 
of course, redundant). That is, loop and x "summations" are done simultaneously. This 
construction, i.e. generation of loops at x = 1 and stochastic summation over x by random 
translations, automatically reproduces (statistically) the correct loop multiplicities in the 
hopping parameter expansion, by which a certain loop of length L appears L times, unless 
it represents r windings of a loop of length L', L = r L' , in which case it only appears 
V times []l2f , |l4| . To improve the loop ensemble we can implement random rotations 
of the loops (before translations). Since the loop evaluation is vectorized over rotations 
and translations, these are reasonably cheap. In the following we work in D = 3 and 
D = 4 and use an {u>(x^ = 1)} ensemble of up to 600000 (250000) loops on lattices of 
N x x N y x N z (xNt) with = 12, 16, 24 and 32, and periodic boundary conditions. 

For a free theory (Ui = 1 above) the contribution of the RW-loops of length L is pro- 
portional to their number, (which counts the different loops with their multiplicities). 
From (111) we have: 

n L ~ T-°l 2 oc L- D ' 2 (8) 
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Figure 1: Number of closed loops of length L, ni generated by random walk. We show the 
frequencies i>l ^) vs L in D = 3 (diamonds) and D = 4 dimensions (crosses). The full 
lines represent the corresponding continuum results L~ D I 2 . 

for large L - this well known relation can be seen by choosing A = 0, m = in (|l]) and 
rescaling r — ► tT, x^ — > x^y/T. In Fig. [l] we show the ul distribution in D = 3 and 
4 for our working ensembles of loops. We observe the excellent agreement with (||) for 
large L and the expected systematic deviations for small L. In the following we shall 
use the actually realized ni (instead of L~ D I 2 ) to normalize the loop contributions when 
necessary, to improve on the systematic errors. 

Notice that the number of loops at large L is much smaller than at small L. If large 
loops are physically important this fact may lead to large statistical errors. One improves 
on this by using a larger number of trials at high L, producing in this way more loops 
there. Then we must renormalize their contributions according to the number of trials. 
The loop frequencies are thus: 

vl= jz t^§7TT? = L-°/ 2 , Ml = nr. of trials at L (9) 
(2tt/D) d / 2 M l L>> i 

(the factor 2 appears because on the cubic lattices we only have even loop lengths). 

3 Special configurations 

For a first test of our approach we use special gauge field configurations: a constant 
magnetic field in full or half space and a constant topological charge density. 
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3.1 Homogeneous fields 

For homogeneous field configurations the continuum path integral can be analytically eval- 
uated to give the well known Euler-Heisenberg-Schwinger Lagrangian (here in Euclidean 
space-time): 



C(c,b) 



1 



(Att) d / 2 



T sinh(cT) sinh(6T) 



(10) 



where 



E 2 + B 2 - \/(E 2 + B 2 ) 2 - 4(EB) 2 



E 2 + B 2 + J(E 2 + B 2 ) 2 - A(EB) 2 
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Of course before performing the T integration one has to renormalize, i.e. to subtract the 
integrand at T = for D = 3 and to subtract also its T 2 dependence for D = 4 



1 9 
1 - -x 2 



sinh x sinh x 6 
The constant magnetic field can be realized on the lattice as 



U n ,i 



exp 



n„ 



1. 



,Ny=Ny, 



U 



71.2 



exp —2irin Q 



, n x = 1, . . . , N x = N x , 



all other links being 1. This provides a field 

b x = 2nk x /N x , b y = 2irk y /N y , b = b x + b y , 

where k x , k y are integers, as required by the periodic boundary conditions. 
In Figs. [2] and || (left hand plot) we show the quantity: 



(12) 

(13) 
(14) 

(15) 



" L 



( 



F(L,b) 



(16) 



i=l 



for various lattices. From (10,11]) (with B z = b, all other fields zero) this should be 

F c (L;b) 
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u 



bL 
2D' 



sinh -u 

For small fields the points fall well together and reproduce the continuum result (11 



(17) 



They also scale well with the lattice size - which, by (13|-15|), fixes the roughness of the 
potential, hence indirectly the lattice spacing. For large fields deviations are seen on small 
lattices, which should be explained by lattice artifacts and signal the departure from the 
continuum. Notice that the data in the figures are plotted against bL, therefore for large 
b we see the small L contributions. 
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Figure 2: Contributions of loops of given length, F(L, b) vs x = bL, in a homogeneous 
magnetic field b for various b in units of it/8 on a 32 3 lattice in D = 3 (left plot) and in 
units of ir / 6 on a 24 4 lattice in D = 4 (right plot). 



We also consider SU(2) configurations of the form: 

U n,l = ex P (^ nia 3 n yJ^ , n y = l,...,N y , 
exp (-2mazn t j^\ , n t = 1, . . . , N t , 



U. 



n,3 



all other links being 1. We take iV„ = N and 

k y = k t = k; b=2nk/N 



(18) 
(19) 

(20) 



where k is an integer. Such a configuration has E = B = (0, 0, b) and hence topological 
charge density 

_ 2^__ti 2 _ 
q ~ N 2 ~ 8tt 2 
homogeneously distributed over the lattice. We calculate 



(21) 



1 UL ( \ 
nL i=i \i e rf J 



(22) 



shown in Fig. |3|, right hand plot. From (|T(],11) this should be 

■ 2 bL 



F c (L,b) 



u 



sinh 2 u ' 



u 



AD' 



(23) 



Again we see good agreement and scaling both with the field and with the lattice size. 
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Figure 3: Contributions of loops of given length, F(L, b) vs x = bL for various lattices N 
in D = 4: in a homogeneous magnetic field b (b in units of ' Air/N) (left plot) and in a field 



of constant topological charge density q 



, b in units of2n/N (right plot). 



3.2 Magnetic field in a half space 

Let us consider now a magnetic field in D = 3 in a half space (x-direction) pointing into 
the z-direction. Take N x = N x /2 in (|i"4| ) and m = 0. We calculate 



£(b, x) 



1 1 nL 

^2-u L F(L,b,x); F(L,b,x) = — ^Tr 



\ 



n ui 

yi G 4(x) j 



(24) 



where vl are the correct loop frequencies (||). Here the loops in Co l L {\) are obtained from 
the original ensemble co l L (l) by approximate centering at the point Xn = 1 (and are hence 
confined in a radius ~ yL around this point). Then the &l(x) are obtained from the 
£b % L {l) by translation with x lattice units in the x— direction. We denote 



1 



C(L 1 ,L 2 ,b,m,x) = ^2 jV L F(L,b,m,x) 



L=Lj 



From the original T— integral (|l|), as discretized in @, we should have: 

Ai A2 



b ' b 



,b, [xVb, £,Vb + 



,T 



6-5/(A 1 ,A 2 ,a*,0- 



(25) 



(26) 



where xq (£ = 0) is the separation between the regions with and without magnetic field. 
We illustrate in Fig. |I] the test of this scaling relation. We find out that it is rather well 
fulfilled if we take xq = + 0.5 - notice, that according to (14,15) the field extends with 
half strength in the region 4^-, + 1 before vanishing. 



(27) 



On physical grounds, see also pi] , we expect 

/(0,oc,M)oce- M ^ 
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Figure 4: Scaling propagator f(X\, A2, /x, £) eg. ^dj ) for /i = 0, various b (in units of tt/16). 
Left plot: f vs £ for various summation intervals Ai — A2- Right plot: f vs X\ for /i = 0, 
A2 = 180 and £ = —5 (non-vanishing field region). 



In the actual calculation we have < Ai 2 < 00 and discretization artifacts, which bring 
systematic errors on M. On Fig. ||, right hand plot, we illustrate the dependence of / on 
Ai; M itself, however, appears less sensitive - see Fig. ||, left hand plot, and Fig. ||. In 
Fig. [5| we show the "effective mass" 



Af Bff (Ai,A 2 ,M,0 = b~ 1/2 In^r^- p ^ (28) 



We observe reasonably good scaling, at least for not too large b and m (where small 
loops dominate and therefore the discretization artifacts are enhanced). Notice that the 
m > results continue nicely to the m = ones, indicating that we do not need to 
take extrapolations but can directly work at m = 0. We read off from these figures 
M(0) = 1.7 =t 0.2. The error is estimated from the sensitivity to the other parameters. 
Since we expect statistical errors to be smaller than this we did not attempt a jack knife 



analysis. Note that our result for M(0) is about a factor 2 larger than that of ref. [11 
This disagreement is beyond the estimated error. 

To understand the roots of the systematic error we observe that in order to perform 
the T integral, we need reliable estimates over the full T range. Consider the simpler case 
of a constant field in the full space: For small T below T = (or, respectively, T = in 
the spirit of a Riemann sum) corresponding to our smallest path with L = 4 the integrand 
factor (bT/ smhbT — 1) can be safely approximated by —^(bT) 2 for bT < 1. For large 
paths with L > Luax = WD/b corresponding to bT > 5 the &T/sinh&T can be neglected 
against -1 in the integrand, i.e. we only need paths in the outer field below Luax- The 
integrand above Tuax = L<Max/2D is easily evaluated (as well as the corresponding L sum 
producing a £ function). Both estimates indicate that the integrand at small and large L 
can be approximated by extrapolation. 
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Figure 5: Effective mass M c g(Ai, A2, //, e#. /or various summation intervals \\ — A2 
and various b (in units ofn/W): vs £ for \i = (7e/f p/o£j and vs [i 2 , £ ~ 0.5 (right plot). 

4 Discussion and Conclusions 

We have introduced the Random Walk Worldline method as a statistical summation of the 
loop expansions for (logarithms of) determinants, and have given examples for effective 
actions of rather simple gauge fields. A more exciting application would be to lattice 
discretized nontrivial field configurations like instantons, sphalerons, bounce solutions. 

We had a scalar particle in the loop but it is very simple to change to fermions and 
gauge bosons, just adding a spin term in the worldline Lagrangian containing the outer 
field |J. In the case of a (nonabelian) gauge field background spin 1/2 corresponds to a 
term ia^F^, spin 1 to 2i F^ (Feynman gauge). These terms add knots to the Wilson 
loop and can be discretized as J7-plaquettes U^U^U^U^ to be joined to the path u>l in 
all possible ways with a a^ v or 2-factor at the conjunction for spin 1/2 or spin 1, respec- 
tively. E.g. the Euler-Heisenberg formula for fermions is thus obtained from the bosonic 
case very easily. There is an overall Dirac and color trace for fermions and a Lorentz and 
color trace for gauge bosons in the loop. The latter of course require C-connections in the 
adjoint representation. The spin term can be further modified by introducing Grassmann 
variables. This is very useful in analytical calculations, in particular if worldline super- 



symmetry is used [17], ||. Also scalar /pseudoscalar and axial vector couplings to spin 1/2 
fermions can be treated that way. 

Alternatively we can stay with the direct expansion of the fermionic (logarithm of the) 
determinant and perform the Dirac traces the way it is done for the hopping parameter 



expansion [14|. 



Even more interesting is the application of the method in the case of a fluctuating gauge 
field background in 3- and 4-dimensional QCD and 3-dimensional thermal electroweak 
theory. This leads to genuine nonperturbative phenomena. 

In the general case the quenched matter free energy is given by @ where we use the 
ensemble of loops generated by random walk. Varying (5 we tune the lattice spacing. The 
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size of the lattice should be chosen such as to accommodate the loops which contribute 
significantly to 

Of course one can discuss strongly interacting gauge theories like QCD fully with 
lattice methods. Still it is a relevant question how the elementary fermion/gauge field 
propagators (here in a loop) are deformed in presence of such backgrounds. The answer 
might allow a semianalytical treatment of some aspects. 

In this context the discussion of chiral symmetry breaking seems to be most promising. 
For our simple case of a constant 5-field this is easy to explain in an analytical treatment: 
calculate e.g. in d = 3 with the Banks-Casher formula for the fermion in the loop 



— d m , N , . f°° dT 1 BT m 2 T 



e 



m=0 



8mB f°° , , r 2 B 

dxe~ x = — (29) 



(47r) 3 / 2 m Jo 2-7T 

(this is well known [[Tgl ). A discretization in case of a fluctuating gauge field background 
appears to be very promising p|. One also may like to develop some modified perturbation 



theory with an IR cutoff like in the stochastic vacuum model |2C] and thus discuss a 



"magnetic gauge boson mass" [||, perhaps with a tachyonic component, in this context. 
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